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X-ray emission from high-redshift miniquasars: 
self-regulating the population of massive black holes 
through global warming 
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ABSTRACT 

Observations of high-redshift quasars at z > 6 imply that supermassive black holes 
(SMBHs) with masses M > 10 9 M Q were in place less than 1 Gyr after the Big 
Bang. If these SMBHs assembled from "seed" BHs left behind by the first stars, then 
they must have accreted gas at close to the Eddington limit during a large fraction 
(<; 50%) of the time. A generic problem with this scenario, however, is that the mass 
density in M ~ 10 6 M Q SMBHs at z ~ 6 already exceeds the locally observed SMBH 
mass density by several orders of magnitude; in order to avoid this overproduction, 
BH seed formation and growth must become significantly less efficient in less massive 
protogalaxies through some form of feedback, while proceeding unabated in the most 
massive galaxies that formed first. Using Monte-Carlo realizations of the merger and 
growth history of BHs, we show that X-rays from the earliest accreting BHs can provide 
such a feedback mechanism, on a global scale. Our calculations paint a self-consistent 
picture of black-hole-made climate change, in which the first miniquasars — among 
them the ancestors of the z ~ 6 quasar SMBHs — globally warm the intergalactic 
medium and suppress the formation and growth of subsequent generations of BHs. 
We present two specific models with global miniquasar feedback that provide excellent 
agreement with recent estimates of the z — 6 SMBH mass function. For each of these 
models, we estimate the rate of BH mergers at z > 6 that could be detected by the 
proposed gravitational- wave observatory eLISA/NGO. 

Key words: black hole physics - cosmology: theory - galaxies: formation - quasars: 
general - gravitational waves 



1 INTRODUCTION 

The discovery of bright quasars at redshift s z > 6 i n the 
Sloan Digital Sky Survey (see the revi ew by Fanll2006r). th e 
Canada-France High-z Quasar Survey (|Willott et al.ll2010al ) . 
and the current redshift recor d-holder at z = 7.08 in the 
UKIDSS |Mortlock et al.ll201ll ) indicates that supermassive 
black holes (SMBHs) as massive as several 1O 9 M0 were al- 
ready in place when the Universe was less than 1 Gyr old. 

The mechanism by which these early massive BHs 
formed and g rew remains poorly understood (see, e.g., 
iHaimanl |2012| for a recent comprehensive review). One 
class of explanations is the rapid collapse of primordial 
gas into a 10 4 — 10 6 Mq black hole, eithe r directly (e.g. 
lHaehnelt fc Reesl Il993l : lLoeb fc Rasid 119941 1. by accreting 



onto a pre—existing smaller seed BH (e.g. IVolonteri fc Reesl 
2005), or by goi ng through the interim st ate of a very 
massive star (e.g. iBond. Arne tt & Carr 1984 ), a rapidly ac- 
cretin g massive "proto-star" ([Begelman. Volonteri fe Reesl 
120061 : iHosokawa. Omukai fc Yorkei l2012h~ or a dense 
stellar cluster JOmukai. Schneider fc Haimanl 120081 ; 
iDevecchi fc Volonteril [2009). These models rely on 
rapid gas contraction in deep potential wells (halos 



temperatu res 
Loebl [20031 



T ■ 



10 4 K; lOh fc Haimanl 



with virial 

| 2002t iBromm fc Loebl 120031: lLodato fc Nataraianl 1200 
Spaans fc Silkll2006l: IWise fc Abelll2008l : iRegan fc Haehnel 



20091 : IShang. Bryan fc Haimanll201of i" and require the gas to 



avoid fragmenting early on by cooling via H2 and/or metals. 
It is unclear whether this special configuration — warm, 
dense, me tal-free gas in relatively massive halos — is realized 
in nature (|Diikstra et al . 2008: ISchleicher. Spaans fc Gloverl 
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Petri. Ferrara fc Salvaterrall2012l ; llnavoshi fc Qmukail 
Agarwal et alj|2012l) . 



An alternative possibility, which we explore further 
in this paper, is that metal-free stars, with masses 
10 — 100 Mp), form at redshifts as high as z > 



25 jAbel, Bryan fc Norrnanll2002l; lBromm 1 Coppi fc Larson! 



120021 ; lYoshida, Omukai fc Hernquistl ^POct). leave b ehind 
remnant BHs with similar masses "(|Heeer et al.ll2003h , and 
subsequently grow b y mergers and by mas s accretion near 
the Eddington limit jHaiman fc Loeb|[200ll ; lHaehnelt| [2003; 
Bromlcv, Somcrvillc & Fabian 2004; Yoo & Miralda-Escude 
20041; lTaniguchill2004l;IShapirdf2005l ; [volonteri fc Reesll2006l ; 
Tanaka fc H aiman 2009), hereafter TH09). The primary the- 



oretical uncertainty of this scenario is whether the seed BHs 
can sustain such high accretion rates. 

To reach masses of > 10 9 Mq by z ~ 7, the stellar- 
mass BHs must grow near the Eddington limit without sig- 
nificant interruption ( with a mean duty cycle > 0.5 at ra- 
diative efficiency 0.07; |TH09)Q However, the accretion rates 
of BHs are limited by their gaseous environments: the rates 
can be significantly sub-Eddington if the seed mass is light, 
or if the BH is embedded in a low-density and/or high- 
temperature medium, as may often be the case at high 
redshifts (see § 12 . 1 . 2 p . Another limitation is radiative feed- 
back: recent work has shown that the accretion is episodic, 
with time-averaged Eddington duty cycle of at most ~ 1/3 
i|Milosavlievic et al. 2009; Par k fc Ricottill2012j : see also ear- 
lier work by ICiotti fc Ostrikerl l200ll for radiative feedback 
forcing accretion to be episodic with a low duty cycle at 
lower redshifts). This local feedback could be mitigated if 
the flow is non-spherical (i.e., in a disk) and the radia- 
tion escapes vertically without compromising the equato- 
rial fuel supply. Finally, a third type of limitation could be 
the lack of continuous fu e l su pply on larger scales. Inter- 
estingly, |Pi Matteo et al.l l|2012T ) have found that, at least 
for > 10 5 Mq BHs residing in massive ( > 10 9 Mq) galax- 
ies at lower redshifts 6 < z < 12, filamentary accretion of 
cool gas into the host galaxy may deliver gas to the central 
regions of massive galaxies near the rates corresponding to 
the BH's Eddington limit. The large Eddingt on ratios (~ 1; 
Willott et aH |2010b) and duty cycles ( > 0.5; IShankar etail 



20101 ) inferred from observations of the quasar SMBHs at 
z ~ 6 appear to be consistent with a sustained and rapid 
BH growth. 

If the BHs can indeed maintain such high accre- 
tion rates, then they can grow into > 1O 9 M SMBHs 
by z > 6. However, another, less appreciated problem 
then arises: these optimistic assumptions inevitably lead to 
a severe over-production for the space density of lower- 
mass (10 5_7 Mq) BHs. The global comoving mass den- 
sity of SMBHs in galactic nuclei in such models can ex- 
ceed the locally observed mass density by several or- 
ders of magnitude, even if l ess than 0.1% of star- forming 
halos form a seed BH ([Bromley. Somerville fc Fabianl 
120041 ; iTanaka fc Haimanl :2009). The overproduction can be 
avoided if the growth of lower-mass nuclear BHs is sup- 
pressed at late times, e.g. by imposing an early M — a rela- 



1 Alternatively, the accretion could occur intermittently in brief 
episodes of super-Eddington accretion. 



tion or scaling the BH growth with the host galaxy merger 
history. 

The purpose of the present paper is to consider an alter- 
n ative so lution to this problem. In the SMBH growth models 
of lTHOQl the z ~ 6 SMBHs originate from seeds born in the 
first minihalos at redshifts z > 25. On the other hand, the 
lower-mass (M & 1O 7 M ) BHs arise primarily from seeds 
born inside minihalos collapsing later, at z ~ 15. In prin- 
ciple, therefore, the overproduction of the > 10 5_7 Mq BHs 
can be avoided, provided that the seed BHs are either un- 
able to form, or unable to grow, in the vast majority of 
minihalos at z ^ 20. A natural reason for this could be 
photoionization heating of the intergalactic medium (IGM) 
by the earliest accreting seed BHs themselves. The X-rays 
emitted from these holes have a long mean free path and es- 
tablish an early X-r ay background, which can pre-heat and 
pre-ionize the IGM fohl|200ll: IVenkatesan, Gfroux fc Shunl 
120011 : iRicotti fc Ostrikerl 120041 ; iMadau et all l2004h . Once 
this heating sufficiently elevates the IGM gas temperature, 
the collapse of gas into low-mass halos will be suppressed 
throughout the Universe 

We examine whether X-rays produced from the earliest 
BHs themselves may provide sufficient heating to avoid the 
overproduction of the 10 5_7 Mq BHs at z ~ 6. Our study 
extends that of TH09, by self-consistently treating global 
X-ray heating and the corresponding rise in the halo mass 
scale for seed formation and BH growth. 

This paper is organized as follows. In Sj2] we describe 
our modeling the co-evolution of BHs ( i|2.1|l and the ther- 
mal and ionization state of the IGM ( <I2.2p . Results of the 
coupled evolution of BH growth and IGM heating are pre- 
sented in Sj3] In i|4] we discuss earlier studies on high-redshift 
miniquasar feedback, and offer general conclusions on viable 
models of the 2 > 6 quasar SMBHs. Finally, we summarize 
our main conclusions in iJU 



2 METHODS 

We couple a merger-tree BH assembly model with radiative- 
transfer calculations of the global heating and ionization of 
the IGM. Our semi-analytic method self-consistently models 
the intertwined evolution of the entire population of accret- 
ing nuclear BHs in halos with masses M > 3 x 10 4 Mq in the 
redshift range 6 < z < 45, with a statistical representation 
equivalent to » 5 comoving Gpc 3 . 



2.1 Black hole growth 

Follo wing ear lier work (e.g.. IVolonteri. Haardt fc Madaul 



2003: 



2004 



2006 



.Bromley. Somerville fc Fabiarj 
Yoo fc Miralda-Escudel |2004| ; IVolonteri fc Reed 
TH09), we compute the hierarchical growth of dark matter 
(DM) halos by means of Monte Carlo merger trees in the 

2 Another possibility is metal-enrichment: once the IGM is pol- 
luted, massive PopIII stars may stop forming, and remnant BHs 
will become much less common. However, this is less attractive, 
because the metals will be much more localized around pre- 
existing sources (at least much more than X-ray photons), and 
cannot produce a global feed back (though clustering can m ake 
metal- feedback more effective; [Kramer, Haiman & Oh 20Q(|). 
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extended Press-Schechter f ormalism (lLacev fc Cold 1 19931 ; 
ISomerville fe KolattJ Il999l ; IZhang. Fakhouri fc Mai l2008h . 
and couple it with a semi-analytical model to follow the 
growth and dynamics of BHs. We employ a ACDM cos- 
mology with the parameter values h — 0.704, Oa = 0.728 , 
fi m = 0.272, tt b = 0.045 and erg = 0.81 Jjarosik et al.ll201lT l. 
We start at z — 6 and model the assembly history of the full 
halo mass function above Mhaio > 10 8 Mq at this redshift, 
following the method of TH09. The reader is referred to 
that paper for a detailed description of the Monte-Carlo 
algorithm; here we summarize the m ost imp ortant features 
and highlight the improvements over ITH09I . 

The rarest, most massive halos are simulated individu- 
ally, which determines the cosmological volume represented 
by our suite of simulations. The population of less mas- 
sive halos is divided into logarithmic mass bins of width 
A log(Mh a io/M ) = 0.5; within each bin, a sufficiently large 
number (~ 10 2 " 3 ) of unique halos are simulated to obtain 
a statistical sample for low-mass halos. The BHs in these 
low- mass halos are then counted multiple times ("cloned") 
to represent the same large volume as for the most massive 
halos. Our most massive bin consists of three halos with 
Afhaio > 10 12 ' 7 Mq, from which we infer that our suite of 
merger trees represents a statistical sample equivalent to a 
comoving volume of ~ 5 Gpc 3 . Our results are based on a 
total of « 10 4 merger trees in each model. 

The most significant change relative to |TH09| is that 
we have improved our halo mass resolution to corre- 
spond to a virial temperature of T v i r = 400K. Above 
this va lue, the gas can collapse through efficient cooling 
by H? jHaiman. Thoul fc Loebl IT9961 ; iTegmark et "all 1 19971 ; 
iMachacek. Bryan fc Abelll200l[) . We follow the merger his- 
tory of each halo down to this limit. As a result, our merger 
trees also extend to higher redshift (as large as z max ~ 45), 
and to lower masses (as low as 3 x 10 4 M )0 

At each redshift, the comoving density pbh(z) of BHs 
residing in (proto-)galactic nuclei changes due to a combina- 
tion of three effects: accretion onto the holes already present 
at redshift z, creation of new seeds in small halos, and ejec- 
tions due to gravitational recoil; 



PBH.net (z) = p\cc{z) + Pcrcat (z) ~ peje(z) 



(1) 



BH accretion deposits high-energy photons into the IGM, 
heating it and affecting the conditions for subsequent BH 
growth. Below, we describe our model implementation of 
each of these components. 



2.1.1 Seed black hole formation 

Starting from the highest redshift in each merger tree, we 
follow all of the branches of the tree towards lower redshifts. 
We place seed BHs in halos when they first reach a virial 
temperature of 400K. This corresponds to a mass threshold 



MtZ d) > 9-1 x 10 4 



-3/2 f \ + , 



21 



-3/2 



Mr: 



(2) 



3 The relative streaming velocities (1 — 10 km s -1 ) be- 
tween baryons and dark matter can increase the mass thresh- 
old for halo virialization and delay PopIII star forma tion 
llTseliakhovich. Barkana fc Hiratal l20l.lt IGreif et all l2011at Li, 
Tanaka and Haiman, in preparation.) 



where p is the mean molecular weight. 

The initial mass function (IMF) of the first stars 
remains highly uncertain. Simulations that include the 
effects of radiative feedback from the accreting pro- 
tostar had suggested that the maximum mas s is ~ 
320 M (|Omukai fc Pallal l200ll ; lOhkubo et all |2009h . 
with recent studies yielding values as low as 30 M 
ilStacv. Greif fc Bromrnl 120121 ) to 43 M© (|Hosokawa et all 
l201ll ) Recent simulations have also suggested that hydrody- 
namical turbulence in star-forming halos may place typical 
masses of PopIII stars significantly lower than previously 
believed, perhaps M* > 10 - 50 M ; that their IMF was 
less steep than PopI stars, and that they may have formed 
in small groups or cl usters (TTurk. Abel fc O'Sheal 120091; 
Stacy. Greif fc Bromrnl |20ld: IGreif et al.ll2011bl; iTurk et all 



20121 ; I Greif et al.ll2012r ). For our purposes, what matters is 



only that there is a reasonable chance for forming at least 
one massive PopIII star in a minihalo that can leave behind 
a seed BH. For concreteness, we chose an IMF with stellar 
masses in the ra nge 20 M < M* < 320 M , with a Salpeter 
ijSalpeterl [l955h power-law slope dn/dlogM, oc M . 
This choice is conservative in that most of the seed BHs 
have low masses - but we emphasize that it is still possible 
that most, or even all, PopIII stars had still lower masses 
and left no seed BHs. 

For simplicity, we also assume that at most one mas- 
sive star forms per halo, because the halo is metal-polluted 
once it undergoes a star-formation episode. (This assump- 
tion is conservative: massive stars outside the pair instability 
SN range are thoug ht to collapse dire ctly into a BH with- 
out ejecting metals: [Heger et all 12003 .) Pristine halos that 
merge with such halos are also considered to be polluted. 
We do not allow seed BH formation in metal-polluted halos; 
we consider such halos to be unable to form massive stars 
whose BH remnan ts have masses of more than several M 
l|Heger et al.ll2003h . 

The seed BH masses are modeled using fitting formula e 
to the simulation results of lZhang. Wooslev fc Hegerl (2008), 
who calculated the remnant masses of metal-free stars un- 
dergoing core-collapse explosions while accounting for the 
mass of the ejecta that falls back onto the remnant. Specif- 
ically, we prescribe 



Mbh = 



■(M.-20 M ) + 2M ifM*<45M 
|(M,-20 M ) 



if M* > 45 M ® 



and assume that stars in the pair-instability mass range 
140 M < M* < 260 M leave no BH remnant. 

For reference, we not e that the above differs from the 
earlier treatment of lTH09l . in which all of the BH seeds were 
assumed to have the same mass of 100 Mq. In the present 
model, the mean PopIII stellar mass is (M t ) ~ 49 M , and 
the mean seed BH mass is (Mbh) ~ 13 M . 



2.1.2 Black hole accretion 

As our fiducial accretion model, we take a simple scenario 
in which all seed BHs grow continuously (at least in a time- 
averaged sense) at a fraction /Edd of the Eddington accretion 
rate, 



M = /EddA/Edd; 



(4) 
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Here 
take 



Me 
e 

to 



= 2.2 x 10" 



-Mbh Myr oc M BH . 



(5) 



the radiative effic 



lcncy, 



which we 



be 0.07 fcf. iMerloni fc Heinz! 120081 : 



IShankar. Weinberg fc Miralda-Escudel 120091 1. With this 
prescription all BHs grow exponentially with an e-folding 
timescale of Myr. We also consider a more realistic 

prescription in which the accretion rate is further limited 
by the Bondi rate M Bon di, i-e. 



M = min(/ Edd M Edd , M Bo ndi) 



Mb 



4ttG 2 A/Ih^ ocMIh- 



(6) 



(7) 



Here p gas and c s are the gas density and sound 
speed, respectively. In general, in the cold and dense 
cores of the cuspy gas de nsity profiles of minihalos, in 
which gas cools efficiently (jAbel. Bryan fc Normanl 120021 ; 
iBromm. Coppi fc Larsor]|2002h . the Bondi rate exce eds the 
Eddington r ate, and so does not impose a limitation IjTurnerl 
1991; TH09). On the other hand, the growing seed BHs can 
find themselves in a much lower-density medium for sev- 
eral reasons. First, the progenitor star could have blown 
most of the gas out of the host halo, so that the seed 
BH i s born in a low-density gas (Whalen. Abel fc Normanl 
|2004 I Alvarez. Wise fc Abell l2009h . Second, although the 
photo-heated gas could cool (via Compton scattering on 
the CMB) and fall back into the halo, it will retain ex- 
cess entropy and may not be able to form H2 and cool 
again, e specially if a strong Lyman- Werner background i s 
present (|Haiman. Rees fc Loebl Il997l : lOh fc Haimar]|2003h . 
In this case, the gas will be in a pres sure-supported config- 
urati o n at much lo wer density (e.g. IShapiro. Iliev fc Ragal 
Il999l ; IRicottil liOQgh H Third, even as the halo undergoes 
mergers with other minihalos, the gas in those partner halos 
may have been unable to cool, as well, owing to the presence 
of a Lyman- Werner background. Fourth, the seed BH can be 
kicked out via gravitational recoil from its parent halo, and 
spend a prolon ged period wandering in the lower-density 
halo outskirts (|Blecha fc Loebl l20od : ITH09I ; iGuedes et at] 
I2009T ). 

Rather than attempting to model the above effects self- 
consistently, we simply adopt a fiducial Bondi rate, scaled to 
the central gas density of the ho st halo. In an NFW DM halo 
(N avarro, Frenk fc Whitel 19961 ), the gas density profile has a 
central c ore that is proportional to the background IGM gas 
density (|Madau. Ferrara fc Reed [200 ll ). p g aT°' oc Pga? M ' oc 



(1 + z) 3 . Following IRicottil (2009), we assume an isothermal 
central gas equation of state and adopt the functional form 



(core) 



<5core — 



Pga: 



(IGM) 



(1-0.2/3/3) (c/0.2) 3 
'■ (c/0.2) 3 ^-°- 2 P/ 3 ) -0.2/3/3' 



(8) 



Above, Ag a a ~ 5 is the mean gas overdensity of the halo; c is 
the halo concentration parameter; /3 = cT / /(c) is defined in 
terms of the parameter F = T v - 1t /Tigm and /(c) = ln(l+c) — 
c/(l + c). An important feature of this model is that the core 
density depends on the IGM temperature (this dependence 

4 Note that progenitor stars with masses above 40 Mq are ex- 
pected to collapse directly into a BH and not eject any metals 
l|Heger et alJl200a ). eliminating the option of cooling via heavier 
elements. 



arises through the external pressure from the IGM near the 
halo's virial radius). In the limit of large /3 oc T v i r /TiGM, 
Score ~ 600c 3 . For the concentration parameter, we adopt 



4.5 



/ M h 



V 10 6 Mg 



-1/8 



l+z 
21 



(9) 



where the normalization and the depe ndences on Mh a io and 
z ar e based on the s i mulat ion results of lBullock et al. (|200ll ) 
and IWechsler et ail (|2002D . Several studies have found that 
the redshift dependence of the concentration parameter 
weakens at high redshifts z>4, an d that c has a charac- 
teristic minimum value Cuin ~ 3 — 5 ijZhao et al.ll2003l 120091 ; 
iKlypin. Truiillo-Gomez fc Primackll201ll ). We impose a con- 
servative lower limit c > 3 to equation [9] this still allows for 
the Bondi accretion rate of our BHs to be many orders of 
magnitude below the Eddington rate. 

To summarize, the core gas density is a function of (i) 
the mean IGM gas density (or redshift), (ii) the halo concen- 
tration parameter c, and (iii) the parameter T = T v i r /TiGM. 
Given the large uncertainties in modeling the relevant gas 
density profiles, this Bondi accretion model can be more 
generally understood in terms of the core overdensity 5 COIC . 
The critical value of <5 C orc at which Maondi = /EddAfEdd is 
given by 



2 x 10" 



Mr 



50 Mr. 



T ■ 

1 vir 

400K 



3/2 



l + z 
21 



(10) 



2.1.3 Black hole mergers and gravitational recoil 

Following the merger of two DM halos that each contain a 
central BH, we assume that the BHs promptly form a bi- 
nary and coalesce , aided by the gas su rro unding them (e.g. 



Escala et al.1l2004l ; ICuadra et aTli2009l ; see iKulkarni fc Loebl 



20121 for the possible consequences on the high-redshift 



SMBH population if merger timescales for SMBHs are longer 
than those for their host halos). A recoil velocity is assigned 
to the merged binary, wit h a random value dr awn from the 
fitting formulae given by lLousto et~aH |2010l) . assuming a 
uniform distribution for the spin magnitudes (in the range 
< a < 0.93) and angles with respect to the binary's or- 
bit al plane (0 < 9 < n/6). The la t ter ch oice i s motivated 
by iBpgdanovic, Reynolds fc Milleil l|2007t) and iDotti et all 
(2009), who have argued that coplanar circumbinary disks 
can cause moderate spin alignment. 

Radial trajectories of the kicked BHs are computed as in 

iTHOgl . including the effects of dynamical friction and assum- 
ing an NFW DM potential superimposed with an isother- 
mal gaseous profile (p ga s oc r~ 2 ). The latter is taken to be 
an isothermal cusp in models with an exponential accretion 
prescription (eq.[5]), and the gas profile with a core (eq. [8jl for 
models with the Bondi-limited prescription (eq. [7]). Instead 
of computing the radial trajectories of recoil explicitly for 
each merger event, we compare the recoil velocity against a 
retention velocity threshold v rct (z, Mhaio, Mbh) which is in- 
terpolated from a precompiled table. We consider a recoiling 
BH to be ejected if it does not return to the halo center on 
a radial orbit within a Hubble time. 

In general, we have found that, somewhat counter- 
intuitively, the final SMBH mass function at z > 6 is not 
highly sensitive to the precise spin prescription l|TH09h . or 
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to the central mass distribution of the host galaxy. This is 
because the recoil velocity is a strongly increasing function 
of the mass ratio M2/M1 < 1, and the effects of gravita- 
tional recoil may be understood rather simply in terms of 
critical mass ratios. Nearly-equal mass BH binaries should 
almost always be ejected, barring an ultra-massive host halo, 
strong gas drag from an extremely dense central gas pro- 
file, o r highly speci fic, recoil-minimizing BH spin orienta- 
tions (Haiman 2004). Conversely, mergers with mass ratios 
Si 1 : 20 result in kicks of only Si 20 — 30 km s _1 regardless 
of BH spins, and cannot eject BHs from deep inside their 
halos. Fully-grown SMBHs are ejected only very rarely, as 
only the most extreme mass ratios and spins can produce 
kicks large enough to remove them from their massive host 
halos. Therefore, the spins and the host gravitational po- 
tential can influence only the fates of BHs that merge with 
mass ratios ~ 1 : 10 in low-mass halos early on. 

In our models, the ancestors of z ~ 6 SMBHs tend to be 
BHs that were sufficiently more massive than their merger 
partners, and whose mergers occurred in sufficiently massive 
host halos. The probability that particularly massive BHs 
survive gravitational recoil events increases with time, as 
they become successively more massive with respect to their 
contemporaries an d the gravitational potential of their host 
halos grow deeper l|Volonteri fc Reesll2006l ; ITH091) . 

2.2 IGM heating and radiative feedback 

2.2.1 (Mini-) quasar emission spectra 

We model all accreting black holes as light-bulb (mini-) 
quasars with Eddington luminosities that are "on" for a frac- 
tion Af/AfBdd of the time, and "off" the rest of the time. We 
assume that a fraction / cor = 0.1 of the total energy is emit- 
ted by a hot X-ray corona, modeled as a power-law spectrum 
above 1 keV with a photon index T — 2.0. The rest of the 
emission i s modeled as a multicolor S hakura-Sunyaev accre- 
tion disk (|Shakura fc Sunvaevlll973h , with M = MEdd and 
viscosity parameter a = 0.1. We model the spectrum as a 
one-zone graybo dy photosphere IjB laes 2004]) using the fit- 
ting formulae of lTanaka fc Menoul (|2010l ) for the graybody 
source function; the disk emission is scaled down by a factor 
1 — /cor at all frequencies to accommodate the coronal emis- 
sion. All miniquasar photons with energies above 13.6 eV are 
assum ed to escape the host halo l|Whalen. Abel fc Norman] 
l2004h . 

Figure \T\ shows two examples of our model spectra, 
for miniquasars with M B h = 100 M Q and 10 5 M Q . For 
the entire mass range of interest, miniquasars have signif- 
icant emission above 1 keV (analogous to high-mass X- 
ray binaries and ultraluminous X-ray sources). Our quali- 
tative findings are robust as long as miniquasars are abun- 
dant in the early Univer se and are prolific producers of ~ 
IkeV X-ray emissi on (cf. IVenkatesan. Giroux fc ShullllioOll : 
iMadau et aHl2004l ; IrlkTotti fc Ostrikeril2004l ). 



2.2.2 IGM heating and cooling 

We treat the IGM as a primordial H+He gas with a mass 
fraction Yke ~ 0.24 in He. Initially, it cools adiabatically as 

2 / -, , \ 2 



Tigm(z) = 2.73(1 +zt) 



l + zt 



K « 9 



1 + z 
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Figure 1. Model spectra of accreting BHs, assumed to be light 
bulbs radiating at their Eddington luminosity. Top: Mbh = 
10 2 M Q . Bottom: M B h = 10 5 M Q . The total spectrum (solid 
black) consists of an o-viscosity accretion disk (dashed blue) and 
an additional power-law corona at > IkeV (dotted red). 



where z t » 140 is the redshi ft at which th e IGM and CMB 
temperatures decouple (e.g., Peebles 1993). 

We follow the temperature and ionization state of the 
IGM from high to low redshift, together with the evolving 
BH population, using small timesteps with Az ~ 0.1. (At 
z > 30, where miniquasar heating is negligible, we use Az = 
1.) Across each timestep, the high-energy photons from the 
accreting BHs heat and ionize the IGM; the conditions in 
the IGM in turn influence the formation and growth of seed 
BHs. The ionizing photon background is built up via the 
cumulative miniquasar emission, with the comoving specific 
luminosity density 



t v (z) 



V 



(12) 



where V — 5 Gpc 3 is the comoving volume of our merger- 
tree simulations. To compute the mean background flux at 
reds hift z, we solve the cosmo logical radiative transfer equa- 
tion l|Haardt fc Madaulfl996T ) 



1 

4% 



, t~'\ p -ru{*,*') ( 1 + Z ) 
€v{Z)e (1 + Z') 



(l + z') S dz 



c-t— dz' , (13) 



where 



t„(z, z) 



J z 



. nl (1 + z") c 



dt 

17' 



1 " 

dz 



(14) 



K, (11) 



is the opacity. Above, rig denotes the comoving number den- 
sity of ion species j (HI, Hel, or Hell), o{, is the photoion- 
ization cross section, and v = v" (\ + z)/(l + z"), where v" 
is the emission frequency at redshift z" . 

With the ionization background given by Eq. (|13[) , we 
follow the thermal evolution of the medium by solving, for 
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each time-step dt, the energy conservation equation: 

du d ( 1 \ A nct , . 

= ~P^7 ' ( 15 ) 

dt dt \ p gas J p gas 

where it, p and p gas are the specific internal energy, pressure 
and density of the gas, respectively, and A nct is the net heat- 
ing/cooling rate per unit volume. Heating includes Comp- 
ton and photoionization heating, while cooling includes line 
and continuum cooling (and Compton cooling once the IGM 
temperature rises above Tomb); the first term on the right 
hand side represents cooling due to adiabatic expansion. 
The energy equation is solved coupled with the chemistry 
equations, solving for the abundances of H, H + , He, He + 
and He ++ . The physics of heating and ionization by soft 
X-ray photons has been described in a number of works 
dHaardt fe Madaulll996l: IVenkatesan. Giroux fc Shulllliool 
iFurlanetto fc Stoeverl l2010l ). X-rays propagate to much 
larger distances than UV photons before being absorbed, 
photo-ionizing He and H atoms. These photoionizations pro- 
duce fast photo-electrons, which then partially photoionizc 
the gas through r epeated secondary ion i zation s. We use the 
recent results by IFurlanetto fc Stoeverl (|2010l ) for the frac- 
tion of the energy of each fast photo-electron that is used 
for ionizations and heating. Specifically, we use simple fitting 
functions to the fractions shown in their Figure 4, and we 
neglect the dependence on the photoelectron's energy, which 
is weak for the relevant range considered here (E > IkeV). 



3 RESULTS: EFFECT OF IGM HEATING ON 
BH GROWTH 

3.1 SMBH assembly without feedback 

To illustrate the overproduction problem, we first show in 
Figure [2] the SMBH growth history in two reference models 
without global feedback. The left-hand panels show the re- 
sults of a model assuming exponential accretion (eq. 2| with 
/Edd = 0.55; the right panels show a second reference model 
assuming Bondi-limited accretion (eq. [6| and /Edd = 0.8. In 
both cases, the BHs do not "know" about the temperature 
evolution of the IGM. The core gas densities in the Bondi 
model are computed through eq. [8] as if the IGM cooled 
adiabatically. 

The top panels show the growth of the mass of the most 
massive BH in the simulation (solid black curves), along with 
the comoving mass density of all BHs (blue dotted curves) 
and of only massive BHs with M > 10 5 M Q (blue dashed 
curves). The most massive BH is a particularly massive seed 
with Mbh ~ 100 M©, left behind by a progenitor star with 
M* ~ 260 M©, just above the pair instability window. Be- 
cause the most massive SMBHs in our models are products 
of multiple mergers, their masses scale roughly linearly with 
the maximum allowed seed mass. 

Given that the lo cal SMBH mass den s ity is Pbh.o ~ 

4 x 10 5 Mq Mpc" 3 (jMarconi et al.l |2004| ; IShankar etall 
120041 ') and that > 90% of this mass density appears to 
have accreted at z < 6 during luminous qua sar phases (e.g., 
IShankar. Weinberg fc Miralda-Escudj|2009l ). the mass den- 
sity of massive BHs with M > 10 5 " 6 M at z ~ 6 should be 
at most ~ few x 10 4 Mq Mpc -3 . Both models presented in 
Figure [2] overpredict the global SMBH mass density by sev- 



eral orders of magnitude, due to numerous seed BHs growing 
exponentially near the Eddington rate. 

The failure of the Bondi-limited accretion model to 
avoid this overproduction is not obvious, and is worth ex- 
plaining in some detail. We find that the seeds surviving 
their recoil events and producing the > 10 9 Mq SMBHs at 
z — 6 are those with initial masses of ~ 30 — 100 M©. As 
can be seen from equations [8] and 1101 these seeds, born in 
halos with T v i r > 400K, quickly get out of their Bondi-phase 
and switch to the Eddington-limited regime, provided that 
the central gas density in the halo is sufficiently large (i.e. 
that the halos have virial temperatures sufficiently above the 
IGM temperature threshold). As a result, the most massive 
BHs in the Bondi-limited and the exponential models are 
similar (compare the black curves in the top left vs. top 
right panel in Figure (5]). As equation 1101 shows, the Bondi- 
limit is important for lower-mass seeds, and it also becomes 
more important at lower redshifts, where the characteris- 
tic densities are lower. This is indeed causes a "dip" and 
slows down the overall growth of the BH population by ac- 
cretion in the range 10 Ss z 20. In this range, pBH is domi- 
nated by relatively low-mass 10 3 M ) BHs, which arise 
from relatively low-mass seeds, and which suffer prolonged 
episodes of Bondi-limited accretion. However, as Figure [5] 
shows, this slow-down is insufficient: the total mass density 
in > 10 5 M© BHs by 2 10 is dominated by relatively 
massive and early- forming seeds, which did not spend a sig- 
nificant time in a Bondi-limited stage. 

The bottom panels of Figure [2] show the growth rate 
Pbh of the comoving BH mass density (solid black lines), 
and its three components: creation of new seed BHs (dot- 
ted green lines), accretion onto existing seeds (short-dashed 
blue lines), and ejections due to gravitational recoil (long- 
dashed red lines). The growth of the total BH density can 
be broadly divided into two distinct epochs, being domi- 
nated by the formation of new seeds at early times (z > 20) 
and by accretion at late times [z 15). The seed formation 
history is identical in the two models f £|2.1,l[) . In the Bondi- 
limited case, only BHs with Mbh exceeding a critical value 
Af cr it fsee imp grow exponentially and the rest accrete at a 
much lower rate. This leads to a wider distribution in BH 
masses at late times, which results in fewer BH mergers with 
comparable masses and thus fewer ejections by gravitational 
recoil. In the exponential growth model, BHs tend to have 
similar growth histories; two BHs that have similar masses 
at a given time will also have similar masses at a later time, 
unless one grows much more rapidly via mergers. Thus, pair- 
ings of BHs of comparable masses are common, resulting in 
much higher rates of BH ejection. In both models, and gener- 
ically for assembly models of this type, ejected BHs tend to 
be somewhat lower-mass BHs residing in lower-mass halos 
( H2.1.3p . As seen in the bottom left panel in Figure the 
global mass density of nuclear BHs can temporarily decrease 
during phases in which the ejections of low-mass BHs out- 
weigh the sum of accretion onto high-mass BHs and the seed 
formation rate. 

In both accretion models, the growth rate of the overall 
population pbh becomes unrealistically high by z ~ 90 For 



5 We have checked that despite these high rates, the 2 > 6 
SMBHs in our models never consume more than 10 percent of 
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Figure 2. The evolution of the BH population in two mod- 
els with different accretion prescriptions without global feed- 
back. Left panels: All BHs accrete exponentially at a fraction 
/Edd = 0.55 of the Eddington limit. Right panels: Accretion is 
limited by the local Bondi rate (see text) and by /Edd = 0.8. 
Top panels: The mass of the most massive BH in the simulation 
(solid black curves), and the comoving density of all BHs (dotted 
blue) and of only BHs with M > 10 5 M Q (dashed blue). Bot- 
tom panels: The growth rate pbh °f the global BH mass density 
(solid black), and contributions to it from seed creation (dot- 
ted green), gas accretion (short-dashed blue) and recoil-ejections 
(long-dashed red). 



reference, to avoid overproduction, this rate should have a 
time-averaged value of 

0.1 X pBH,0 



<Pbh(z > 6)) ^ 



At(z > 6) 



60 M Q Mpc" a Myr . (16) 



3.2 SMBH assembly with global IGM warming 

The most conspicuous consequence of IGM heating is the 
inhibition of the gravitational collapse of gas into low-mass 
halos . The cosmological Jeans mass (e.g., iBarkana fc Loeq 
l200ll ) 



Mj(z) = 



( 5kT 1GM y /2 ( 3 y /2 



6.8 x 10' 



4np h 

3 / 2 / Tigm \ 
V 100K J 



(17) 



3/2 



l + z 
21 



-3/2 



Mr: 



will rise, and gas inside halos whose masses are below this 
threshold wil l not collapse to form stars and seed BHs. It has 
been shown |Gnedmll20o"o1 ) that the so-called filtering mass, 
which depends on the temperature history of the IGM, pro- 
vides a better fit to gas fractions in low-mass halos in nu- 
merical simulations than the instantaneous Jeans mass. The 



the baryons in a halo. The most massive SMBHs typically have 
less than one percent of the host halo's baryon mass. 



filtering mass lags the evolution of the Jeans mass, and is 
therefore lower than the Jeans mass, during the period when 
the IGM temperature is rising. However, we note that the 
definition in equation [17] is a factor of 8 lower than the one 
adopted bv lGnedinl (|2000t ). and is, in fact, very close to the 
filtering mass shown in his Figure 3 at the relevant redshifts 
(z > 10)0 We note that the impact of IGM heating may 
extend beyond the fiducial Jeans mass: even in a halo above 
the Jeans mass, gas falls in, at least initially, along a higher 
adiabat from the pre-heated IGM — this could prevent the 
gas from condensing to a sufficiently high enough de nsity 
to activate H2 coo l ing. [Machacek, Bryan fc Abell (|2003T ) and 
iKuhlen fe Madaul ([2005) both found evidence in 3D simula- 
tions that heating the IGM by X-rays can suppress baryonic 
infall into low-mass halos. 

Another consequence of IGM heating may be to 
suppr ess the BH accretion rate in pre-existing ha- 
los. iMachacek. Bryan fc Abell (|2003T ) and lKuhlen fc Madaul 
( 2005) found that an early soft X-ray background can elevate 
the mean gas temperature inside low-mass minihalos with 
masses below the collapse threshold. (However, both studies 
found that the X-ray feedback on cold gas may be posi- 
tive in more massive minihalos; see i]4.1l below.) As the IGM 
heats, the inflow of cool gas into existing halos whose masses 
fall below the collapse threshold would be suppressed. The 
heating and evaporation of small halos could in turn im- 
pact accretion onto BHs residing in more massive halos, 
by reducing the amount of cold and dense gas supplied by 
minor mergers. Furthermore, the accumulating X-ray back- 
ground can heat the gas inside low-mass halos near or above 
their virial temperatures, redu cing the gas density, or pho- 
toevaporating the hal o entirely l|Shapiro. Iliev fc Raga|[20o3 ; 
lAhn fc Shapircll2007n , shutting down Bondi accretion onto 
the nuclear BH. 

There is some circumstantial evidence that the BH ac- 
cretion rate is tied to the mass supply rate onto the parent 
halo. Massive BHs appear to have been much more active 
in the early Universe: most z ~ 6 quasars appear to be 
shining near or slight ly above their Eddington luminosities 
(|Willott et al.l l2010bl ) , c ompared to typical E d dington ra- 
tios o f ~ 0.1 at z < 4 (iKollmeier et al.l 120061 : iKellv et all 
120101 ). IShankar etail |2010l ) concluded that most SMBHs 
are ac tive as quasars at z ~ 6, compared to less th an 1% at 
z < 2 jShankar, Weinberg fc Miralda-Escudell2009l ). A plau- 
sible physical reason behind such hyperactivity is the rapid 
mass supply to t heir host halos via m ergers with other halos 
and by accretion; lAngulo et alj (|2012l ) observed that numer- 
ically simulated host halos of z ~ 6 quasars doubled in mass 
in the preceding 100 Myr. 

Motivated by these considerations, we simulate two 
models of SMBH assembly that include negative feedback 
due to IGM heating. In both models, the cosmological Jeans 
mass rises due to IGM heating by miniquasars; new seed BHs 
are only formed in chemically pristine halos with masses 
above the Jeans mass. In the first ( "exponential" ) model, 
we assume that BHs grow exponentially with /Edd = 0.55 



6 iNaoz fc IB arkana ( 2007|) showed that the gas density and pres- 
sure gradients in the initial conditions after recombination de- 
creases the filtering mass. However, this is below our minimum 
halo mass imposed by cooling, and does not affect our results. 
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as in one of the reference models, but that this growth is 
enabled by the supply of cold gas from the merger-driven 
growth of the host halo. Thus, we halt BH accretion one 
e- folding ti me (» 62 Myr , comparable to a typical quasar 
lifetime; e.g. lMartinill2004l ) after the most recent merger of 
their host halo with another halo. We additionally require 
that the merger partner halo exceeds the Jeans mass (so 
that it contains cold gas). In the second ( "Bondi" ) model, 
we assume that BH accretion is instead limited by the Bondi 
rate in the halo nucleus, which we model as in the second 
reference model. We again take / E dd = 0.8 and p gas = p m i n . 
However, the central gas density of the halo depends on the 
ratio r = Tvir/TiGM (eq. [8} - as the miniquasar emission 
heats the IGM, the core densities in low-mass halos decrease 
and the Bondi accretion rates of their BHs are suppressed. 



3.2.1 Global IGM warming: results 

Our results in both the Exponential and the Bondi models 
are shown in Figure [3] As in Figure [2] the left panels show 
the exponential accretion model, while the right panels show 
the Bondi-limited model. The top panels show the growth 
of the most massive BH in each simulation, along with the 
global BH mass density Pbh, while the bottom panels show 
Pbh and its components due to seed formation, gas accretion 
and recoil-ejections. The color and line style schemes are 
identical to those in Figure [2] For ease of comparison with 
the reference models, the vertical scales of the graphs are 
the same as in Figure [2j despite much lower values of pbh 
in the models with feedback. 

In the top panels of Figure [3] we show the z = 6 BH 
mass function produced by our models, without (thin lines) 
and with feedback (thick lines) . As with the previous two fig- 
ures, left panels show the exponential-accretion models and 
the right panels show the Bondi-limited models. Feedback 
clearly depresses the late growth of Mbh ~ 10 5-7 M© BHs, 
while leaving the most massive objects nearly unaffected. 
Note that the mass function is very steep in the models 
without feedback. This is because most seeds form within a 
relatively brief cosmological epoch, and a large number of 
BHs are allowed to grow at the identical exponential rate. 

In the bottom panels of Figure [3] we estimate the lu- 
minosity function for the models with feedback that avoid 
SMBH overproduction. Recall that in order to treat the 
X-ray emission from accreting BHs, we have approximated 
them as Eddington-luminosity lightbulbs with a duty cycle 
M/MEdd- Using this prescription to translate the mass func- 
tion into a quasar luminosity function, however, will tend 
to overestimate quasar luminosities and underestimate their 
counts (duty cycles). We therefore also consider the oppo- 
site extreme estimate for the luminosity function, in which 
BHs shine at a bolometric luminosity (M/MEdd)-^Edd and 
have a duty cycle of 1. The true luminosity function should 
lie within the bounds set by these two formulations, rep- 
resented in Figure [4] by the shaded grey region. For refer- 
ence, we plot the observational determination of the z ~ 6 
quasar luminosity fu nction (corrected for obscuration) of 
Willott et al. 1 i|2010bl ; blue dotted lines), as well as the fittin 
formula suggested bv lHopkins. Richards fc HernauistJ | 200 
blue dashed lines). 

The luminosity functions predicted by our models, in 
particular the "exponential" model, may still be too steep 
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Figure 3. Same as Figure[2] but global IGM heating by accreting 
BHs suppress BH formation and growth at late times. Left pan- 
els: BHs accrete exponentially as in the left panels of Figure [2] 
but accretion is halted if host halos do not continue merging with 
massive halos containing cold, collapsed gas. Right panels: BH 
accretion is limited by the local Bondi rate as in the right pan- 
els of Figure [2] the rising IGM temperature lowers the halo gas 
density (eq. [S}. In both models, halos below the Jeans mass are 
disqualified from forming new seed BHs. 



below L ~ 10 46,5 erg s -1 compared to the best-fit luminos- 
ity function of lWillott et al.l (|2010bh . However, the faint end 
of the z ss 6 quasar l uminosity function remains highly un- 
certain. The IWillott et al.l fit includes a single binned data 
point below 2 x 10 46 erg s" 1 , and assumes that the frac- 
tion of o bscured quasars (which is larger for low-luminosity 
objects; iLawrencg 1 199ll ) is the same as the one found in 
the local Universe. The latter point is motivated by the fact 
that the obsc ured fraction doe s not appear to evolve between 
z = and 2 l|Ueda et al.ll2003h . but there has been no direct 
determination of the obscured fraction at z ~ 6. 

As Figures [3] and U show, both of these models avoid 
the overproduction of the lower-mass (M ~ 10 5 " 7 ) BHs, 
while still producing M > 10 9 M Q by z Ri 6.4. Note 
that the exponential growth model, in which the accretion 
rate is capped at half of the Eddington limit, falls some- 
what short of explaining the z = 7.08 quasar S MBH with 
M ~ 2 x 10 9 Mq observed by IMortlock etHI i|201lf ) de- 
spite produci n g a go od fit to the z « 6 mass function of 
IWillott et ail l|2010rj ). At z = 7.08, in a simulation vol- 
ume roughly 1/3 of the ~ 15 comoving Gpc covered by 
the UKIDSS, the most massive SMBH in this model is a 
factor ft* 4 short. The z — 7.08 SMBH can be reproduced 
with only a slightly higher accretion rate, /Edd = 0.6; how- 
ever, this overshoots the global mass function and pbh by 
roughly an additional order of magnitude. This illustrates 
the fact that this SMBH is yet more difficult to accommo- 
date in simple SMBH assembly models. One obvious pos- 
sibility that the IMortlock et al.l SMBH is a mild outlier in 
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Figure 4. Top panels: The 2 = 6 BH mass function in our 
SMBH assembly models. The thin histograms show the models 
without feedback, severely overproducing the global SMBH mass 
density; the thick histograms show the models with feedback. The 
error bars demarcate the Poisson errors from our simulated sam- 
ple. Bottom panels: The 2 = 6 quasar luminosity function in our 
models with feedback, shaded in grey between the upper and 
lower limits for the luminosity and duty cycles of individual BHs 
(see text). For reference, we show the obscuration-corrected lumi- 
nosity function found by Willott et al. (2010; blue dashed curve) 
and the best-fit model for 2 = 6 of Hopkins et al. (2007; blue 
dotted curve). As with Figs. [2] and [3] the left panels show expo- 
nential accretion models with /Edd = 0.55, and the right panels 
show models where accretion is capped by the Bondi rate and by 
/Edd = 0.8. 



its growth history. While /Edd = 0.55 is a suitable mean 
value in the context of this model, realistically BHs have a 
distribution of Eddington ratios and not a single universal 
value. 

In summary, these results of our illustrative models sug- 
gest a generic, self-consistent scenario of black-hole-made 
"climate change", in which the earliest and most massive 
accreting BHs heat the IGM and suppress the formation 
and growth of subsequent generations of BHs. As in the 
case of global warming on Earth, the negative effects are 
felt by the next generations: the BHs originally responsible 
for this feedback are largely unaffected by it, because they 
reside in the most massive halos that merge frequently with 
other massive halos and have high central gas densities to 
facilitate BH growth. 



3.2.2 IGM global warming and reionization 

We show the thermal and ionization histories of the IGM 
in Figure [5] As with the previous figures, the exponential 
accretion models are plotted on the left panels, while the 
Bondi-limited models are shown on the right. The mod- 
els with and without feedback are plotted with thick and 



thin curves, respectively. Note that in the models without 
feedback, we still allow miniquasars to heat and ionize the 
IGM, but BH formation and accretion are assumed to be un- 
affected by this cosmic climate change. In particular, seed 
BHs continue to form at M^ d) (T vir = 400K), and BH ac- 
cretion rates proceed universally at /Edd in the exponential 
model, while proceeding as if the IGM were cooling adia- 
batically in the Bondi model. In the models with feedback, 
seed formation and accretion are suppressed by the warm- 
ing of the IGM, as described above in §3.2\ The top pair of 
panels in Figure [5] show the IGM temperature Tigm [z) for 
each model. The middle pair show the corresponding cosmo- 



logical Jeans mass; the seeding mass threshold M l 



(seed) 
halo 



due 



to H2 cooling is plotted in blue dashed lines alongside the 
Jeans mass. The bottom panels show the electron fraction 
x e = (nmi + riHeii + 2nHeiii)/(nH + 2n Hc ). 

Together, Figures [3] and [5] shed light on the coupling 
of the BH accretion history with the IGM temperature his- 
tory. Initially, the evolution of Tigm (2) and Mj(z) follow 
a power-law decay due to adiabatic cooling. The buildup 
of the X-ray background heats the IGM, causing a reduc- 
tion in the seed BH formation rate once Mj becomes larger 
than M^°^ . The suppression of seed formation occurs at a 
unique characteristic temperature, which may be estimated 
analytically as follows. The Jeans mass scales with the IGM 
temperature and density in the same way that the halo mass 
scales with the virial temperature and characteristic internal 
gas density (both require the sound-crossing time to equal 
the dynamical time): M oc (T 3 /p) 1//2 . Accounting for the 
gas overdensity inside collapsed halos and differences in the 
conventional normalization coefficients of order unity, this 



N3/2 



i.e. M. 



AL 



(seed) 
halo 



leads to Mj/M ha i « 60(T lGM /T vi 
when Tigm ~ 26K. The rising Jeans mass quenches seed BH 
formation at z ~ 18 for both the exponential model and the 
Bondi-limited model. At these redshifts, the total nuclear 
BH mass density is already pbh ~ 10 4 MqMpc -3 ; however, 
most (> 90%) of this density still consists of the original 
seed BHs and there are no holes above 10 4 M©. 

We may also understand the energetics of the global 
IGM warming, as follows. The increase in the IGM thermal 
energy density, (3/2)nbfcATiGM — where n D is the baryon 
number density of the IGM and ATigm is the difference be- 
tween the actual IGM temperature and the value expected 
from adiabatic cooling (eq. 1 1 If) — should scale with the total 
energy density emitted by accreting BHs. The latter is equal 
to eAp acc c 2 , where Ap aC c is the global mass density accreted 
by BHs. Taking into account that only a fraction /heat of the 
accretion power goes toward heating the IGM, we obtain 



ATigm « 10 2 /i 



heat 



Apace 

M Mpc -3 



(^) (JL) K . 

V0.07/ V1.2/ 



(18) 



We find that initially, at z > 20, only / hea t ~ 10~ 4 - 10~ 3 
of the (mini-) quasar energy goes into heating the IGM, but 
that this increases to 10 -2 at z ^ 8. The redshift dependence 
of the heating efficiency is due to the fact that not all of 
the hard photons are immediately absorbed by the IGM 
due to their long mean free paths, which results in a delay 
between emission and absorption. Further, only 1 — 10% of 
the spectrum is emitted at E ~ 1 keV, with more massive 
BHs emitting softer spectra. Some of the energy is also lost 
to ionization and to the CMB through Compton cooling. 
From the accretion history of Pbh, we can also con- 
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elude that merely turning off the seed BH formation rate 
is insufficient to prevent overproduction. This is because 
X-ray heating is initially inefficient, with only a fraction 
/heat ~ 10 -4 of the BH rest-mass energy at z ~ 20 going 
into elevating the gas temperature. The global BH density 
must therefore grow to a value as large as 10 4 M Q Mpc -3 
for global X-ray feedback to halt seed BH formation (at 
Tigm>30K). At the redshifts where this occurs (z ~ 20), 
the most massive BHs are at most ~ 10 3 Mg. If all BHs 
were to grow exponentially from this point onward, by the 
time the most massive BHs reach 10 9 M Q , the entire popu- 
lation must grow by similar factors of ~ 10 6 . Even allowing 
for some of the BH growth to occur via mergers, extreme 
overproduction of Pbh is inevitable. To prevent overproduc- 
tion, either another mechanism must regulate (much more 
severely) the formation of seed BHs, or BH accretion must 
be (preferentially) extremely inefficient inside in low-mass 
halos at later times. 

In both models with feedback, the BH accretion rate is 
indeed reduced preferentially inside low-mass halos. In the 
exponential accretion model, BH growth is reduced in host 
halos that do not continue merging with other halos that 
exceed the Jeans mass. (Note that this suppression kicks in 
at z » 25, before seed BH formation rates are affected.) 
In the Bondi-limited model, IGM heating lowers the central 
halo gas densities, and the corresponding Bondi rate, in the 
lower-mass halos. 

Unsurprisingly, the same growth models without feed- 
back that overproduce the massive BH population also ap- 
pear to overpredict its contribution to reionization. Current 
observational constraint s indicate that reionizat i on wa s still 
ongoing at z > 7 (e.g. iFan. Carilli fc Keatind |2006) , and 
that the optical depth of the IGM to Thomson scattering is 
t cs = 0.088 ± 0.014, only about half of which is accounted 
for by full y ionized gas betwee n < z < 6 (CMB mea- 
surements; |Komaisu~etj2[201l]). Our "exponential" model 
without feedback reionizes the Universe by z r = 8.4 and pre- 
dicts an IGM Thomson optical depth of r cs ~ 0.09, while our 
"Bondi" model without feedback completes reionization by 
z T = 9.1 and predicts r cs ~ 0.11. These models reionize the 
Universe earlier than the empirical constraints, despite ne- 
glecting contributions to reionization by the stellar UV radi- 
ation accompanying the mini-quasars BHs. The models with 
feedback, on the other hand, fall well short of the empirical 
constraints (neither completes reionization by z — 6) and al- 
low room for stars to be the prim ary sources of reionization. 
iDiikstra. Haiman &: Loe"bl (|2004h showed that even scenar- 
ios where the contribution to reionization by miniquasars 
is ~ 10 — 20% are consistent with the observed unresolved 
soft X-ray background near 1 — 2 keV. As our models with 
miniquasar feedback have lower contributions to reionization 
(and also because our miniquasar spectra are softer than the 
one used by those authors), we conclude that they are well 
within th e empirical constrai nts on the X-ray background 
(see, e.g.. IDiikstra et ai1l2012l . for a recent study). 



3.2.3 BH merger rates 

In Figure [Bl we show the number of BH mergers at z > 
6 in the binary mass window 10 4 Mq < M(l + z) < 
10 7 Mq, the estimated sensitivity window for major mergers 



Exponential accretion Bondi limited 

40 30 20 15 10 8 6 ^0 30 20 15 10 8 6 
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7, 

Figure 5. Top panels: The IGM temperature as a function 
of redshift. Middle: The cosmological Jeans mass (black solid 
curves) and our fiducial halo mass threshold for seed BH forma- 
tion (dashed blue curves). Bottom: The IGM electron fraction. 
Thin curves represent models without feedback, whose BH mass 
evolution is given in Fig. [2] thick curves represent the models with 
feedback shown in Fig. \3\ As in Figs. [2]-|4] the exponential accre- 
tion models are shown on the left and the Bondi-limited models 
on the right. 

for the proposed gravitatio nal-wave detector eLISA /NGO 
|Amaro-Seoane et al.ll2012T ). The merger rates are shown in 
our two feedback models, in the previous figures. Both mod- 
els predict that w 30 major mergers (binary mass ratios 
0.1 < q < 1) per year per unit redshift may be detectable at 
6 < z < 10, with app roxima tely twice as many events with 
mass ratios q > 0.01. TH09 found that in a class of simpli- 
fied models that prevent overproduction by requiring that 
the z > 6 quasar SMBHs form from extremely rare seeds, 
the detection rate may be zero over the mission lifetime. It 
is encouraging that in the new models presented here, the 
seeds are systematically common and the detection rates are 
nonzero (see explanation for this difference below). 

The majority of mergers occur with mass rat ios 0.01 < 
q < 0.1 (see also iHollev-Bockelmann et al.l [20101 ). Mergers 
with more nearly equal masses are relatively rare, because (i) 
feedback tends to increase the mass discrepancy of BH pairs 
by preferentially suppressing the growth of lower-mass BHs, 
(ii) BHs that survive recoil events tend to become increas- 
ingly more massive with respect to their contemporaries, as 
explained above, and (iii) halos with very different masses 
are considered unable to coales ce in ou r merger tree, due to 
the long merger timescales fsee iTHolSI for details). 

3.2.4 Effects of the Lyman- Werner background 

We have also run models that included a global Lyman- 
Werner (LW) background produced by the progenitor stars 
of the seed BHs, which can raise the minimum halo mass 
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Figure 6. BH merger rates for the two assembly models with 
feedback in the approximate detection window of the proposed de- 
tector eLISA/NGO, binary mass 10 4 M < M(l + z) < 10 7 M . 
The line styles demarcate the binary BH mass ratio q < 1; in 
order of increasing thickness: any q > 0, q > 0.01, and q > 0.1. 



threshold for subsequent seed formation. This background 
and the associated mass threshold can be calculated as fol- 
lows. Approximately 10 4 photons with energies E w 12.9 eV 
are produced per baryon of a progenitor star. Through 
our Salpeter mass function and our progenitor-to-seed mass 
relation (eq. [3]), the comoving mass density of progeni- 
tors is related to that of the seeds as pp op m ~ 3.91 x 
/° Pcrcat(z) dz. The fraction of the photon energy that 
contribu tes to the background, after radiative transfer, is 
^ 1/ 10 l|Haiman, Rees fc Loeblll997l ; iHaiman. Abel fc Reed 
2000). This gives a characteristic early LW background in- 
tensity of 



J LW « 6 x 10 



PPopIII 



10 4 Mq Mpc" 3 
x erg s " cm 



21 



Hz" 



(19) 



Using the fitting formula f rom the simulations of 
IMachacek, Bryan fc AbeH (|200ll ), this leads to a minimum 
halo mass scale for subsequent seed formation: 



M LW « 2 x 10 b 



PPopIII 



10 4 M Q Mpc" 



1/3 



l + Z 

21 



-1/2 



Mr: 



(20) 

In our simulations, this mass scale is only a factor of a few 
above the fiducial T v i r = 400K threshold (eq. [2)| during the 
brief interval 20 Ss z 25, before the X-ray global warming 
abruptly turns off seed formation (Fig. [3J . 

The net effect of the LW background was negligible in 
the runs that included both LW and X-ray feedback. As 
such, and given the order-unity errors associated with LW 
radiative transfer ("sawtooth" modulation; see Haiman et 
al. references above) and in the merger tree itself, we have 



shown above only the results with the X-ray feedback to 
best highlight the effects of the latter. 

The above calculation may underestimate the early LW 
background, as it does not account for the possibility of 
PopIII seeds forming in groups alongside the progenitors 
of the BH seeds. At later times, halos with virial temper- 
atures above 10 4 K will efficiently form stars and add to 
the LW background; however, this occurs long after the X- 
rays have turned off seed formation. Even if the LW back- 
ground is much higher than what we have estimated here, 
it still may not suppress BH ac cretion inside previously 
se eded halos, due to self-shielding jDraine fc Bertoldil 1 19961 : 
cf . IWolcott-Green. Haiman fc Brvanll201ll and refs. therein) . 



4 DISCUSSION 

4.1 Comparison to Other Work 

X-ray heating and early reionization by miniquasars has 
been studied by a number of authors. The broader ef- 
fec ts of vari ous high-z X-ray backgrounds was su rveyed 
bvlohl (|200ll ) and lVenkatesan. Giroux fc Shulll (|200ll ), with 
a special emphasis on the temperature and ionization 
structure of the IGM before reionization is complete. 
IMachacek, Bryan fc Abel (|2003D used numerical simulations 
to study the effects of an early soft X-ray background on 
the temperature and density of gas inside minihalos at 
20 < z < 30. At these redshifts, the 1 - 2 keV flux in our 
models is comparable to their parameter choices ex = 1 
and ex = 10, and their results are consistent with what we 
find at the onset of feedback in our models. They report 
that X-ray heating results in a moderate reduction (by a 
factor of a few) of the cold gas fraction in minihalos with 
Afhaio & 10 6 Mq, a mass range comparable to the heating- 
elevated Jeans mass in our models. They also found that 
the mean gas temperature inside less massive halos with 
Afhaio £s 10 5 Mq can be elevated by an order of magnitude, 
consistent with our assumption that an X-ray background 
can act to suppress BH accretion in halos with Mhaio Ss Mj. 

iMadau et all (|2004l ) made a first attempt at specifically 
treating the impact of the first BHs on the global IGM in 
more detail, but did not explicitly includ e feedback by X-ray 
heatin g on the later generations of BHs. iRicotti fc Ostriker] 
i|2004h developed a similar semi-analytical model for early 
"pre-ionization" by seed BHs, with the aim of identifying 
a BH growth scenario which could account for the large 
Thomson-scattering optical depth inferred by the earlier 
CMB measurements. Their study, while including the effect 
of IGM heating on suppressing virialization below the Jeans 
mass, did not investigate the co rresponding modification of 
the re sulting BH mass function. IRicotti. Ostriker fc Gnedinl 
(2005) studied X-ray pre-ionization by means of cosmolog- 
ical simulations. However, their focus was again on the ob- 
servational signatures in the CMB and 21-cm signal, rather 
than on the growth of SMBHs. 

More recently, iDevecchi et al.l (|2012l ) explored the for- 
mation of nuclear star clusters and BH seeds, tracking 
the chemical, radiative and mechanical feedback of stars 
on the baryonic component of the evolving halos. In par- 
ticular, they examined the role of the LW photons in 
suppressing star formation, but did not include the rise 
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of the Jeans mass due to X-ray heating of the IGM by 
the seeds BHs. Feedback effects from the LW photons in 
shaping BH growth in the most massive (10 12 ~ 13 M Q ) ha- 
los at z = 6 have also been r ecently exa mined in detail 
by Petri. Ferrara fc Salvaterral (2012) and Agarwal ct al. 
<|2012h . using merger trees and N-body simulations, respec- 
tively; however, those studies did not examine feedback on 
the lower-mass BHs in smaller halos. In principle, LW radia- 
tion could help alleviate the SMBH overprediction problem, 
by disabling the cooling in low-mass halos at very early times 
|Haiman fc Brvarj|2006l ) - this is similar to the X-ray back- 
grou nd, although likely restricted to lower-mass minihalos 
(e.g., IMachacek, Bryan fc Abeill200lfl . 

Finally, X-rays from miniquasars can also exert positive 
feed back by triggering the formation of molecular hydro- 
gen jHaiman. Abel fc Reed l200ol: IMachacek. Bryan fc Abell 
120031 ; IKuhlen fc Madaul 120051; jjeon et al.l |2012|). For ex- 



ample , IMachacek. Bryan fc Abel ( 20031 '). IKuhlen fc Madaul 
(|2005h and Ijeon et al.l (|2012r 7found that the X-rays can 
increase the amount of cold and dense gas (and thus star 
formation) in the cores of massive minihalos in their simula- 
tions. This runs contrary to the assumption in our "Bondi- 
limited" model that X-rays cause a warm ing of the cores of 
low-mass halosO] However, the results of IKuhlen fc Madaul 
(2005) suggest that at sufficiently high X-ray fluxes, the 
heating outweighs the enhanced H2 cooling and the feedback 
on the cold gas density becomes negative. The value of the 
X-ray flux where this turnover occurs in their simulation — 
in minihalos in close proximity of their miniquasar; ex ~ 
a few at ~ 1 keV in t he no rmalization convention of 
IMachacek. Bryan fc Abel l|2003h — is reached in our models 
at z ~ 25 - 20 when p B H > 10 3 M Mpc -3 , and exceeded 
soon thereafter. This suggests that the X-ray background 
fluxes in our models are large enough to induce negative 
feedback. Furthermore, the E > 1 keV X-rays have mean 
free paths of ~ 1 Gpc. The simulatio n vol umes of ~ 1 Mpc 3 
employed bv lKuhlen fc Madaul i|2005h and ljeon et"aH l|2012l ) 
do not account for the X-rays originating from sources out- 
side the box and thus do n ot address the early bu ildup of 
a global X-ray background. (|Kuhlen fc Madau 2005 turn on 
a single miniquasar at z = 21.) Finally, even if the X-rays 
exert a positive feedback by increasing the H2 cooling, the 
net effect may also be negative in the presence of a strong 
LW field accompanying the X-rays. Whether, and to what 
degree, X-rays can induce positive feedback on structure for- 
mation and early BH growth remains very much an open 
question. 



4.2 General Comments on Viable SMBH Growth 
Scenarios 

The puzzle of the origins of the z w 6 — 7 quasar SMBHs has 
motivated numerous recent theoretical investigations. Here, 



remark on possible SMBH assembly models, attempting to 
generalize the results we obtained above. 



4-2.1 Explaining the monsters 

First, what is necessary to explain the individual SMBHs? 
We find that for the most massive BHs in our simulations 
to grow to M ~ 10 9 Mq by z ~ 7 (i.e., to explain the 
iMortlock et~ai1 (|201lh quasar observation), a seed BH left 
behind by a PopIII star and starting at 100 Mq at z ~ 40 
must have an effective e-folding time of ^ 55 Myr. That is, 
the mean Eddington ratio must be /Bdd ^ 0.6 for radiative 
efficiency e — 0.07. Or, more generally, these two quantities 
must satisfy /Edd(l — e )/ e ^8. As the elapsed time between 
z ~ 40 and z m 7.08 is ss 700 Myr, this lower limit corre- 
sponds to accretion-driven growth by a factor of a few x 10 5 . 
Our models show that mergers can contribute to assembly 
by an additional factor 10 — 100. Note that for the progeni- 
tors of the z ~ 6 quasar SMBHs, this requires that the mean 
accretion be close to , that the radiative efficiency be no 
larger than ~ 0.1, or both. Note that the above requirement 
suggests that the progenitors of the z > 6 quasar SMBHs 
have (averaged over their growth histories) accretion rates 
close to or exceeding the Eddington limit, relatively low ra- 
diative efficiencies e ^ 0.1, or both. Such low radiative ef- 
ficiencies suggest that the BHs cannot be spinning rapidly, 
consistent with accretion dis k models with magnetohydro- 
dynamic turbulence (see, e.g. IShapiro| [2005. and references 
therein). A hypothetical, single BH accreting at the Edding- 
ton limit and e — 0.07 can grow by a factor ~ 10 9 in the 
same redshift interval. 

Another possibility is that the ancestors of the mon- 
ster SMBHs were born with much higher masses. For ex- 
ample, particularly massive BHs with Mbh > 10 4 Mq could 
have formed in halos whose g aseous components have unusu- 
ally l o w angular momentum (fKoushiappas. Bullock fc Dekell 
120041 ; lLodato fc Nataraianl l2006h ; in halos that are heav- 
ily irradiated by massive neighbors |Diikstra et al.l [2008); 
in environments with high ambient magnetic fields 
jSethi. Haiman fc Pandevll2010l ); or in rare instances of ion- 
izing shocks in dense halo cores (|lnavoshi fc Omukaillioil ). 

A common minimum requirement for both classes of 
seed models is that a ~ 10 J Mq BH be in place by z « 10; 
such BHs can grow at the Eddington rate to ~ 1 9 Mq by 
z ~ 7. The recent work bv lDi Matteo et al.l (|2012T l suggests 
that once such a BH is in place, then cold gas accretion by 
the rapidly growing host halo could help deliver the gas, at 
least on large scales, to sustain its near-Eddington growth. 
Any seed model that satisfies this condition can, in princi- 
ple, explain the individual quasar SMBH masses observed 
to date. The above requirement also suggest a fundamental 
degeneracy for assembly models of high-redshift SMBHs in 
the 6 Si z 10 range. 



7 If the bottleneck on BH accretion is not the central cold gas den- 
sity but rather the external supply of g as, then the X-ray-induced 
suppr e ssion of baryonic infall fo und by IMachacek. Bryan fc Abell 
j2003h , [Kuhlen fc Madaul j2005l l and others should lead to a neg- 
ative feedback on BH growth, as assumed in our "exponential 
growth" model. 



4-2.2 Modeling the overall population 

In light of the above points, the greater challenge ar- 
guably lies in distinguishing between the various assem- 
bly scenarios. Probing the non-degenerate parameter space 
(Mbh < 10 s Mq,z>10) through direct electromagnetic 
observation will be difficult. This is about the detection 
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sensitivity of the James Webb Space Telescope, assuming 
Eddington luminosity and a ~ 10 s s integration l|Haimanl 
120121 ). If IGM preheating suppresses star formation in ha- 
los with Mhaio ^ 10 9 Mq, this may lead to an observ- 
able faint-en d cutoff in the high-r edshift galaxy luminos- 
ity function l|Barkana fc L ocb 200 0]) and a drop in the rate 
of hi gh-redshift supernovae IjMesinger. Johnson fe Haimanl 
120061 ); both effects could be observed by the JWST and 
provide circumstantial evidence for negative feedback on 
low-mass halos. However, the most obvious candidate for 
distinguishing models is a gravitational-wave detector such 
as eLISA/~NGO, which could directly constrain BH merger 
rates, masses and abundances at high-redshift. 

Previous studies had suggested that to 
avoid overproduction, seed BHs must be rare 
jBromlev, Somerville fc Fabianl 12004 lTH09l ). In prin- 
ciple, the occupation fraction of SMBHs in galaxies can 
approach unity by low reds hift even if nuclear BHs were 
very rare at early times ([Menou. Haiman fc Naravananl 
l200ll ; iLippai, Frei fc Haimanl |2009h . If seeds are indeed 
rare, then it is conceivable that the BH merger rate at high 
redshift z > 6) would be so low t hat no detectable events 
will be expected by eLISA/NGO |TH09l ). A rare seed BH 
population would also have a smaller impact via global 
feedback. Therefore, one way to indirectly confirm such 
models is by ruling out alternative scenarios, with common 
seed BHs, through null detections. 

The alternative possibility proposed here is that seeds 
are systematically much more common, but that BHs grow 
much less rapidly in low-mass halos at late times. This 
could be because AGN feedback or other local processes self- 
regulate the BH mass with respect to the host galaxy proper- 
ties, i.e. processes simil ar to those resulting in the locally ob - 
served M-a relation |TH09l . rVolonteri fc Nataraian| [2009) . 
In this paper, we have shown that similar effects can be pro- 
duced instead by global feedback; the hard photons from the 
early accreting BHs can heat the IGM and regulate subse- 
quent BH formation and growth. The same BHs that go on 
to assemble the z ~ 6 SMBHs can slow the growth of the 
low end of the BH mass function and avoid overproduction. 
One important consequence of the common-seed scenarios, 
regardless of whether the mechanism regulating BH growth 
is local or global, is that the mergers of the SMBH an- 
cestors and their intermediate-mass contemporaries should 
be frequent enough to be detectable by an observatory 
such as e£/&4/N GO fcf. ISesana. Volonteri fc Haardtll2007l . 
IMicic et al.ll2007l l. 

Additionally, if the seeds are common, then the X-rays 
emitted in the course of their mass growth could generically 
influence cosmological structures. This may be inevitable, 
regardless of whether the first nuclear BHs had masses of 
10 M or 10 5 M Q ; en route to becoming 10 7 " 9 M SMBHs, 
the vast majority of the mass growth must occur through 
accretion. Indeed. [T51 suggests that the IGM will be heated 
to ~ 10 4 K if the comoving SMBH mass density at z = 6 is 
several percent of the present-day value, and if ~ 0.1% of 
the miniquasar energy output required to build up this mass 
density goes toward heating the IGM. 

Our models predict that miniquasars play a sig- 
nificant role in establishing the cosmological tempera- 
ture history of the IGM, since X-rays from miniquasars 
would heat the early Universe much more effectively 



than UV radiation from PopHI star s l|Madau et al.l 120041 ; 
ICiardi. Salvaterra fc Di Matted l2010T i. Our results suggest 
that PopIII star formation may be affected by miniquasars 
as early as z « 20. Other sources, such as LW radiation 
from PopIII stars and the accreti ng BHs, hard photons from 
non-nuclear stellar-mass BHs JWheeler fc Johnsonl l201lh 



and high-mass X-ray binaries |Mirabel et al. 



2011), or local 



AGN feedback and regulation, may also have contributed to 
this feedback. 



5 CONCLUSIONS 

Observations of high-redshift quasars at z > 6 imply that 
SMBHs with masses M > 10 9 M were in place as early as 
z > 7. In this paper, we considered models in which these 
SMBHs grow from stellar-mass seed BHs forming at ultra- 
high redshifts (i.e. the remnants of first-generation stars at 
z ~ 30 — 40), obeying the Eddington limit on the mass 
accretion rate. Previous work has shown that this is feasible, 
but only if the duty cycle for accretion is of order unity. 

Here we highlighted a relatively under-appreciated 
problem in this class of models: unless the growth of BHs 
in low-mass halos is preferentially and severely suppressed, 
the models overpredict the abundance of 10 5_7 Mq BHs in 
galactic nuclei by several orders of magnitude. 

Using Monte-Carlo realizations of the merger and 
growth history of BHs, we show that the X-rays emitted by 
the earliest accreting BHs can heat the IGM, and suppress 
the formation and growth of subsequent generations of BHs 
in low-mass halos. In this "global warming" scenario, the 
BHs originally responsible for the warming are largely unaf- 
fected by it, because they reside in the most massive halos, 
well above the Jeans mass, and frequently merge with other 
massive halos with cold gas, facilitating BH growth. How- 
ever, the negative effects are felt by the next generations of 
low-mass halos, in which seed formation and BH accretion 
are suppressed. 

We presented specific models with global miniquasar 
feedback that provide excellent agreement with recent es- 
timates of the z = 6 SMBH mass function. These models 
could be constrained through direct observations by JWST , 
and through the detection of tens of BH mergers at z > 6 by 
the proposed gravitational- wave observatory eLISA/NGO. 

We compared our work to previous studies investigating 
the effects of X-ray feedback on early structure formation. 
A primary uncertainty is whether, and to what degree, a 
moderate X-ray background may induce positive feedback 
on the cold gas content of minihalos through enhanced H2 
formation and cooling, and whether this effect outweighs the 
heating of the gas by the X-rays. Finally, a limitation of our 
merger tree approach is that we do not consider local effects 
due to clustering and t o proximity to miniqua sars, which 
may be significant (e.g.. iKuhlen fc Madaull2005l ). While our 
work serves as a proof-of-concept that global X-ray feedback 
from accreting BHs can regulate their growth, more detailed 
studies on this subject are warranted. 
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